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摘 要 :表面 活性 剂 分 子 溶 于 水 自 组 装 形成 的 蠕虫 状 胶 束 溶液 由 于 其 高 表面 活性 和 高 黏 弹性 ,具有 
独特 的 力学 响应 特性 ,呈现 出 与 牛顿 流体 过 异 的 流动 特征 。 采 用 考虑 溶液 中 胶 束 可 逆 断 裂 和 重 构 
效应 的 双 组 分 VCM 本 构 方 程 , 研 究 了 溶液 流 变 学 特性 ,发 现 溶液 拉 伸 应 力 随 应 变 率 呈现 出 拉 伸 硬 
化 - 稀 化 -硬化 的 强 非 线性 特征 。 研 究 了 蠕 叶 状 胶 束 溶液 绕 圆柱 流动 的 特性 ,在 低 雷 诺 数 (Re ~ 
0.01) 下 , 当 维 森 贝 格 数 Wi» 1 时 , 流 场 出 现 较 强 的 不 稳定 性 , 绕 流 尾 迹 区 流 场 呈现 随时 间 波 动 的 特 
征 ,通过 功率 谱 分 析 发 现 流 场 波动 具有 拟 周 期 特征 , 当 胶 束 断 裂 引 起 流体 拉 伸 应 力 对 弹性 响应 的 非 
线性 增强 时 , 流 场 波 动 的 不 规则 性 增强 。 圆 柱 后 缘 附 近 出 现 一 对 周期 性 变化 的 对 称 涡 ,表明 流体 弹 
性 会 促进 流动 分 离 的 发 生 , 流 体 拉 伸 应 力 对 弹性 变化 的 强 非 线 性 响应 是 造成 低 震 诺 数 绕 流 流动 不 
稳定 性 的 根本 原因 ,揭示 了 蠕虫 状 胶东 溶液 在 低 雷 诺 数 高 黏 弹性 时 绕 流 流动 不 稳定 性 的 产生 机 理 。 
关键 词 :表面 活性 剂 ;蠕虫 状 胶 束 溶液 ;VCM 模型 ;流动 不 稳定 性 ; 黏 弹性 流体 
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Abstract: Wormlike micellar fluids formed by dissolving surfactants in water have unique mechanical 
responses due to its high surfactivity and strong viscoelasticity. These properties enable wormlike micellar 
fluids rich and complicated flow behaviors. By considering the breakage of long chains and the reforming of 
short chains,a two species VCM model is introduced in this study to numerically investigate the effect of 


the chain scission on the flow instability. First, the rheological properties are obtained by directly solving 
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the VCM constitutive equation set numerically. The fluid exhibits a strong nonlinear thickening-thinning- 


thickening phenomenon under the uniaxial tension. The flow passing by a confined cylinder is investigated. 


The flow becomes significantly unstable at very low Reynolds numbers (Re ~0.01) when the Weissenberg 


number is larger than 1. The fluid velocity in the downstream wake of the cylinder presents a fluctuation 


with the time. The flow fluctuation is quasi-periodic according to the power spectrum analysis. The flow 


fluctuation becomes more irregular when the nonlinearity of the variation of the tension stress with the 


elasticity becomes stronger. A pair of symmetrical vortexes behind the cylinder changing with time are also 


observed, indicating that the elasticity promotes the flow separation on the cylinder surface. The flow 


instability is due to the strong nonlinear response of the tension stress within the bulk fluid to the elasticity of 


the fluid. 


Key words surfactant ; wormlike micellar solution; VCM model ;flow instability ; viscoelastic fluid 


表面 活性 剂 分 子 的 水 溶液 在 激活 剂 作用 下 自 组 
装 形成 具有 一 定 长 度 (可 达 几 个 微米 ) 能够 自由 弯 
曲 和 拉 伸 的 柔性 聚集 体 ,又 称 蠕 虫 状 胶 束 , 当 水 中 表 
面 活性 剂 浓度 大 于 交 短 浓度 时 ,蠕虫 状 胶 束 交友 缠 
绕 形 成 三 维 的 空间 网 络 结构 ,具有 显著 的 黏 弹 性 和 
丰富 的 流动 现象 ,广泛 应 用 于 提高 油气 采 收 率 '"、 
强化 传 热 中 减 阻 “ 生物 医药 "环境 修复 和 日 
FETU ,近年 来 成 为 石油 .化 工 和 新 材料 的 研究 
dos, 

A ripe IR CAE EL RC Ps D] B E A E s 8 8 
性 ,与 高 分 子 聚 合 物 溶液 中 胶东 的 强 共 价 骨架 连接 
方式 不 同 ,蠕虫 状 胶东 溶液 中 胶 柬 会 发 生 持续 可 逆 
的 断裂 和 重 构 "? ,因此 又 称 “ 活 性 聚合 物 ”, 其 流 变 
性 质 较 一 般 聚合 物 溶液 更 为 复杂 ,如 在 大 变形 的 前 
切 和 拉 伸 下 胶 束 会 加 速 断 多 ”1 。 

目前 对 蠕虫 状 表 面 活性 剂 淤 液 的 研究 主要 集中 
在 其 流 变 学 性 质 , 如 在 稳 态 或 非 稳 态 剪 切 和 拉 伸 作 
用 下 的 应 力 -应 变 率 响应 特性 "中 和 表面 活性 剂 溶 
液 的 化 学 组 成 对 其 自 组 装 形 貌 及 溶液 流 变 性 质 的 影 
响 "71。 在 实际 流动 中 ,流体 经 常 处 于 复杂 几何 和 
变形 的 共同 作用 下 ,呈现 出 丰富 的 流动 现象 ,如 剪 切 
诱导 结构 , 低 雷诺 数 绕 流 过 程 中 被 绕 流 物体 前 后 的 
ADIERA A ,颗粒 在 溶液 中 沉降 速度 的 随 
机 脉动 等 ”2 “2  。 但 是 目前 对 于 上 述 流动 现象 及 
其 机 理 均 认识 不 足 , 对 流动 状态 难以 预测 ,甚至 不 能 
定性 预测 。 

流体 绕 圆柱 流动 由 于 同时 伴随 非 均匀 前 切 以 及 
拉 伸 变形 而 呈现 丰富 的 流动 状态 ,因此 常常 作为 一 
个 基准 问题 研究 流动 问题 ,其 用 于 研究 复杂 流体 时 
比 单纯 均匀 剪 切 和 拉 伸 流动 能 更 多 地 反映 流体 流 变 
学 性 质 玉 2 2 。 对 蠕虫 状 胶 柬 流体 的 圆柱 绕 流 实验 


发 现 ,即使 在 很 低 雷 诺 数 下 , 当 弹 性 较 强 时 , 绕 圆 柱 
流动 过 程 中 产生 显著 的 流动 不 稳定 ”22 ,说 明 流 
体 黏 弹性 会 显著 影响 流体 动力 学 特性 。 

目前 对 于 复杂 黏 弹性 流体 绕 流 数值 模拟 研究 大 
多 基于 Maxwell 流体 ,如 Oldroyd-B 以 及 考虑 剪 切 稀 
化 效应 的 GIESEKUSP?! FENE 本 构 模 型 "3351 ,这 
类 流体 在 描述 聚合 物 溶液 时 适用 性 较 好 ,但 是 对 于 
表面 活性 剂 形成 的 蠕虫 状 胶 束 活性 大 分 子 , 难 以 描 
述 缠 结 胶 束 链 在 剪 切 和 拉 伸 作用 下 的 断裂 作用 , 采 
用 以 上 本 构 方 程 对 蠕虫 状 胶 束 流体 的 很 多 流动 现象 
难以 预测 ,如 难以 解释 实验 中 颗粒 在 蠕虫 状 胶 束 溶 
液 中 沉降 速度 连续 不 稳定 脉动 现象 。VASQUEZ 
45 1 2077 JEF Cates 的 “活性 聚合 物 ” 理论 提出 双 组 
分 VCM 模型 ,考虑 流体 变形 过 程 中 胶 束 的 可 逆 断 裂 
和 重 构 作 用 ,能 够 更 真实 地 描述 蠕虫 状 表 面 活性 剂 
溶液 的 流 变 性 质 ,模型 能 够 更 准确 地 描述 流动 现象 。 

本 研究 拟 采 用 双 组 分 VCM 本 构 模 型 描述 蠕虫 
状 胶 束 溶液 的 流 变 性 质 ,研究 其 绕 圆柱 流动 时 的 流 
动 状态 ,流体 流 变 参数 对 流动 状态 以 及 流动 不 稳定 
性 的 影响 ,以 期 为 认识 上 述 复杂 流动 现象 提供 理论 
解释 ,为 工程 问题 中 复杂 流体 流动 设计 和 优化 提供 
理论 基础 。 


1 数学 模型 和 数值 模拟 方法 


蠕虫 状 胶 束 溶液 为 不 可 压缩 流体 ,假定 流动 为 
绝热 状态 。 流 体 中 除了 溶剂 的 黏 性 应 力 外 ,溶质 中 
胶 束 的 黏 弹性 产生 一 个 附加 的 应 力 项 ,附加 的 应 力 
项 采用 双 组 分 VCM 黏 弹性 本 构 方程 。 

1.1 流动 控制 方程 


连续 方程 
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VE = 0 (1) 
动量 方程 
p(w tu vu) =V( -pI +us(Vu+(Vu) ) +7) 
(2) 


式 中 :po 是 流体 密度 ;u 是 速度 矢量 ;t 是 时 间 sus RE 
剂 黏度 ;z 为 溶质 产生 的 附加 黏 弹性 应 力 , 需 要 引入 
额外 的 方程 进行 封闭 ,此 处 涂 质 的 黏 弹性 采用 双 组 
分 VCM 模型 描述 。 

流体 为 表面 活性 蠕虫 状 胶 束 溶液 ,微观 尺度 的 
胶 束 链 在 不 可 压缩 牛顿 流体 的 溶剂 中 均匀 分 布 ,并 
且 相 互 绰 结 形成 活性 的 空间 网 络 结构 ,网 络 中 的 链 
在 流动 过 程 中 可 以 变形 .断裂 以 及 重 构 ,这 一 效应 通 
过 引入 双 组 分 VCM 本 构 模 型 进行 模 化 。 活 性 网 络 
中 的 弹性 链 产生 的 附加 黏 弹性 应 力 通过 引入 长 链 组 
分 A 和 短 链 组 分 B. 的 数量 密度 和 构象 张 量 求解 , 表 


T = GJ[ (A *2B) - (n, * ny)IEJ (3) 
AP : 6UALBesR REB EE n LA 分 别 为 长 链 组 分 A 
的 数量 密度 和 构象 张 量 ;ns 和 B 分 别 为 短 链 组 分 B 
的 数量 密度 和 构象 张 量 ;7 为 单位 张 量 。 

VCM 模型 假定 在 流体 流动 过 程 中 , 拉 伸 变形 时 
一 个 长 度 为 上 的 链 A 可 以 断裂 形成 两 个 长 度 相 等 的 
长 度 为 1/2 的 短 链 B ,松弛 后 两 个 短 链 B 可 以 重 构 
形成 一 个 长 链 A ,其 断裂 和 重 构 的 速率 由 链 断 裂 和 
重 构 反应 动力 学 考虑 ,组 分 A 和 B 的 数量 密度 输 运 
方程 为 


on 1 > CAR 

Fa +u: Vn, = 2D, Vin, + De a (4) 
ð 2 
Te pu- Yny =2D, Vins - D +24 (5) 


ot Àa Àa 
RP: DA DSZIHIU IK BE A 和 短 链 B 的 扩散 系数 ; 
和 A 为 长 链 A 的 松弛 时 间 ;cx 和 ca 分 别 为 链 A 和 B 的 
断裂 和 重 构 速率 。 

长 链 A 断裂 形成 两 个 短 链 B 的 速率 根据 当地 
应 变 率 、 胶 束 构 象 张 量 和 数量 密度 确定 , 即 

CA = Ciim T RE - (6) 

RP: ch。 为 平衡 时 的 长 链 A 的 断裂 速率 ,六 = 
1 + cy iy 为 速度 梯度 张 量 , 胶 束 断 裂 因子 二 表征 
TUR BE A 断裂 的 非 线 性 效应 大 小 。 两 个 短 链 B 
重 构 形 成 长 链 A 的 速率 假定 为 常数 ,等 于 平衡 态 重 
构 速率 co = cu 。 
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两 种 链 的 构象 张 量 方程 分 别 为 
jA +Á- nl — A.D. V/A = cynyB -cÅ 


(7) 

M Ng 2 
eAB +B- 51- 8Dr V Bz -2ec,nyB +2ec,A 
(8) 


式 中 € =À Àp AREE B 的 松弛 时 间 。 
无 量 纲 参 数 雷诺 数 Re 和 魏 森 贝 格 数 Wi 分 别 定 
义 为 


Re = (9) 
Mo 
AU; 

Wi = (10) 
Àa 

As 7 puc 2 


式 中 :U Hr IE VERE UI Wi 表示 的 是 弹性 
与 猪 隆 的 相对 大 小 ;A_ 为 流体 有 效 松 好 时 间 。 


1.2 数值 模拟 方法 


研究 对 象 为 低 雷 诺 数 (Re 10 7) 时 蠕虫 状 胶东 
溶液 在 高 为 #H 的 通道 内 绕 直 径 为 d 的 圆柱 的 流动 特 
性 ,如 图 1 所 示 。 考 虑 不 同 的 通道 阻塞 比 BR( BR = 
d/H) ,通道 入 口 段 和 出 口 段 的 长 度 分 别 为 40d 和 
60d ,保证 流体 在 绕 流 前 后 达到 充分 发 展 和 松弛 状 
态 。 对 计算 域 采 用 结构 化 网 格 , 对 于 圆柱 表面 速度 
和 应 力 梯 度 较 大 的 区 域 进行 加 密 , 尾 迹 区 流体 由 于 
剪 切 和 拉 伸 作用 较 强 ,对 其 进行 加 密 , 如 图 1 所 示 。 


1 计算 域 及 网 格 
Fig.1 The computational domain and grid 
采用 OpenFOAM 开源 软件 进行 求解 ,守恒 方程 
和 本 构 方程 的 扩散 项 和 梯度 项 均 采 用 中 心 差分 格 
式 ,为 了 提高 迭代 收敛 性 , 所 有 对 流 项 均 采 用 
CUBISTA ( convergent and universally bounded interpo- 
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lation scheme for treatment of advection ) 离散 方法 , 速 
度 压 力 耦 合 采用 SIMPLE 算法 。 

流动 控制 方程 边界 条 件 为 :在 圆柱 和 通道 壁面 
处 ,速度 采用 无 滑 移 边界 u 2u(0,0,0) ,压力 采用 零 
梯度 条 件 Vp =0; 通 道 入 口 处 采用 速度 和信 口 wu = 
u(U;,,0,0) ,压力 采用 零 梯度 条 件 YP =0; 通 道 出 口 
处 采用 压力 出 口 条 件 = Pu ,速度 梯度 为 0。 

对 VOCM 本 构 方程 ,根据 文献 [11 ] ,边界 条 件 
为 :圆柱 和 通道 壁面 处 , 胶 束 构象 张 量 和 胶 束 数量 密 
度 均 采 用 无 通 量 边界 n. Vo 2006 代表 构象 张 量 
分 量 和 胶东 数量 密度 ) ;通道 入 口 处 , 胶 束 数量 密度 
n, 2 l,n, =V2Ca Cus, =1. 93, 构象 张 量 A = 
(100101),B=(n,/20 0 n,/2 0 n,/2) pŠ E 
口 处 , 胶 柬 数量 密度 和 构象 张 量 均 采用 无 通 量 边界 
n: Vozo, 

在 圆柱 壁面 和 整个 计算 域 采 用 不 同 单元 数量 的 
网 格 进行 网 格 无 关 性 验证 ,经 验证 计算 域 具有 
69 000 个 四 边 形 单元 的 网 格 (圆柱 表面 网 格 数 量 为 
360 ,网 格 尺寸 为 0.008 72) 可 以 很 好 地 捕捉 圆柱 壁 
面 附近 及 整个 计算 域内 的 流动 特征 ,进一步 加 密 网 
格 对 计算 结果 影响 可 以 忽略 不 计 。 为 保证 Courant 
数 小 于 0.1 和 计算 收敛 , 取 计 算 时 间 步 长 At =5 x 
10“s, 胶 束 构象 张 量 分 量 和 数量 密度 初始 值 均 取 通 
道 入口 边界 值 。 蠕 虫 状 胶 束 流体 不 仅 在 剪 切 变形 条 
件 下 表现 出 稀 化 ,而 且 在 拉 伸 变形 中 表现 出 拉 伸 硬 
化 - 稀 化 效应 ,因此 选取 参数 必须 能 够 同时 反应 这 两 
种 流 变 特性 ,选取 的 流 变 参数 为 B = wu, = 0.26, 
1. 22.6,2 20.005, n 21.93, 


mi 


2 结果 和 讨论 


2.1 VCM 流体 流 变 性 质 


为 了 研究 VCM 流体 流 变 性 质 ,分 别 考虑 流体 在 
均匀 平面 前 切 流动 和 单 轴 拉 伸 流 动 下 的 应 力 响应 特 
性 。 考 虑 到 VCM 本 构 方程 自身 的 非 线性 ,难以 获得 
解析 解 , 可 通过 给 定 剪 切 和 拉 伸 变形 率 ,然后 直接 数 
值 求解 双 组 分 VCM 本 构 方 程 (3) ~ 8) 获得 流 变 曲 
线 ,结果 分 别 如 图 2 ~ 图 3 所 示 。 

均匀 剪 切 变形 时 流 变 曲线 如 图 2 所 示 。 横 坐标 
Wi =Asy, "y 为 剪 切 速率 ,有 效 松 弛 时 间 As 根据 式 
(11) 求 得 为 0.38 s。 图 2(a) 和 图 2(p) 纵 坐标 分 别 
为 无 量 纲 剪 切 应 力 rw 和 剪 切 黏度 no "TELE HS, 
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VCM 模型 可 以 描述 溶液 的 前 切 稀 化 效应 和 前 切 带 ， 
这 与 Giesekus 本 构 模 型 ” 剪 切 性 质 类 似 ,此 外 , 随 着 
胶 束 断裂 指数 YX = 人 ) 逐 渐 减 小 , 剪 切 带 对 应 的 
剪 切 速率 逐渐 后 移 , 当 X 减 小 至 零 时 , 与 Oldroyd-B 
本 构 模 型 的 剪 切 流 变 曲线 一 致 , 剪 切 稀 化 效应 消失 。 
胶 束 断裂 指数 X 反映 流体 剪 切 应 力 随 弹 性 大 小 变化 
的 非 线 性 强 弱 程度 ,X 越 小 , 剪 切 应 力 随 弹性 越 接近 
于 线性 变化 。 


107 107 10° 10' 10° 10° 10* 
Wi 
(a) 剪 切 应 力 曲 线 


10! 2 3 10* 


(b) 850] 5 BE Hf 4X 


2.8%( 质 量 )OTAC+0.56%( 质 量 )NaSal 


(c) VCM 模 型 预测 值 与 流 变 实验 结果 对 比 


图 2 VOM 流体 在 平面 剪 切 变形 下 的 流 变 曲 线 (网 =Ary ) 
Fig.2 The rheological properties for VCM fluids 


under the planar shear deformation (Wi=Auy ) 
采用 控制 应 力 流 变 仪 ( HAKKE-RS300 ) 对 表面 
活性 剂 溶 液 (2. 8wt% OTAC 和 0.56wt% NaSal) 的 
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2(c) 所 示 。 溶 液 松 弛 时 间 A 20.27 s, 零 剪 切 黏度 
Ho =1.944 Pa . s, 采 用 VCM 本 构 模 型 对 实验 数据 进 
行人 参数 拟 合 , 拟 合 结果 如 图 2(c) 所 示 。 结 果 表 明胶 
RRI AX 为 1 时 稳 态 剪 切 黏度 的 VCM 模型 预 
测 值 跟 实验 结果 吻合 较 好 。 
单 轴 拉 伸 变 形 时 流体 流 变 曲线 如 图 3 所 示 , 横 
坐标 Wig = Asae, 2 为 拉 伸 速率 ,图 3(a) 和 图 3(b) 
纵 坐 标 分 别 为 无 量 纲 拉 伸 应 力 rr 和 拉 伸 黏度 ma。 
可 以 看 出 拉 伸 应 力 随 着 拉 伸 速率 的 增加 呈现 很 强 的 
非 线 性 , 呈现 拉 伸 硬化 - 稀 化 -硬化 的 特征 ,不 同 于 
Giesekus 流体 只 出 现 拉 伸 硬 化 的 情况 ;胶东 断裂 
指数 X 同样 反映 拉 伸 应 力 随 弹 性 变化 的 非 线 性 效 
应 光 越 小 , 拉 伸 应 力 随 弹 性 变化 的 非 线性 越 强 , 随 
TEX YÉCK , 拉 伸 应 力 随 弹 性 变化 非 线 性 逐渐 减弱 。 
因此 , 维 森 贝 格 数 表征 弹性 大 小 ,而 胶 束 断裂 指数 表 
征 流体 应 JI 随 弹 性 变化 的 非 线 性 强 弱 程度 。 


x10* 


X107 


10° 107 10° 10' 10° 
Wi; 
(a) 拉 伸 应 力 曲线 


107 10" 10° 10' 10 
Wi, 
(b) 拉 伸 黏度 曲线 
图 3 VOM 流体 单 轴 拉 伸 变 形 下 的 流 变 曲线 ( Wis =Aure ) 
Fig.3 The rheological properties for VCM fluids under the 


uniaxial tension deformation ( Wi: — A,g& ) 


2.2 VCM 流体 圆柱 绕 流 流动 特性 


通过 对 蠕虫 状 胶 束 浴 液 圆柱 绕 流 的 流 场 特征 研 
究 , 发 现在 低 雷 诺 数 (Re -10 77) Wi » 1 时 尾 迹 区 速 
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度 场 出 现 显著 的 不 稳定 性 ,如 图 4(b) 所 示 , 而 低 雷 
诺 数 下 牛顿 流体 圆柱 绕 流 流 场 呈 稳 态 状态 (图 4a)， 
此 外 采用 其 他 黏 弹性 本 构 方 程 如 Oldroyd-B、 
Giesekus 模型 时 ”5 预测 的 流 场 也 呈 稳 态 流动 ， 
说 明 流 体 弹 性 并 不 必然 导致 流动 的 不 稳定 。 图 4 
(c) X 1 210.45 s 时 圆柱 前 后 对 称 轴 上 无 量 纲 速 度 
分 布 。 可 以 看 出 ,牛顿 流体 绕 流 前 后 速度 分 布 关 于 
圆柱 中 心 对 称 ,而 VCM 流体 绕 圆 柱 流动 前 后 速度 分 
布 在 圆柱 上 游 与 牛顿 流体 类 似 , 而 在 圆柱 后 缘 附 近 
很 小 的 一 个 区 域内 速度 出 现 了 一 个 负 值 , 即 流体 出 
现 与 来 流 相 反 的 负 速 度 ,形成 负 尾 迹 ,这 与 采用 
Giesekus 模型 研究 结果 '”| 类似, 负 尾 迹 是 弹性 和 前 
切 稀 化 共同 作用 形成 。 

为 了 获得 颗粒 下 游 速度 定量 变化 规律 , 取 流 场 
m X(x/d) =1,Y(y/H) =0 点 的 速度 随时 间 的 变 
化 ,如 图 4(d) 所 示 , 可 以 看 出 ,此 处 速度 大 小 随时 间 
的 变化 呈现 近 于 周期 性 的 波动 规律 。 对 得 到 的 速度 
波动 信号 提取 其 功率 谱 。 图 5(a) 和 (b) 所 示 分 别 为 
不 同 现时 速度 变化 时 间 谱 和 功率 谱 , 可 以 看 出 , 当 
Wi =1 时 ,速度 随时 间 并 不 发 生变 化 ,而 当 下 增加 
至 2 时 ,速度 随时 间 呈 现 较 强 的 波动 性 ,其 无 量 纲 速 
度 波动 幅度 在 0.05 ~0.75。 由 图 5(b) 的 功率 谱 可 
以 看 出 ,速度 波动 是 由 数 个 不 同 频率 ,不 同 能 量 的 波 
动 信号 又 加 而 成 ,有 一 个 起 主导 作用 的 主 频 信号 , 主 
频 较 低 , 信 号 频率 越 高 ,幅度 越 小 ,说 明 速 度 在 以 一 
个 主 频 波动 过 程 中 还 有 比较 小 的 信号 的 扰动 ,因此 
绕 流 流 场 速度 的 波动 是 一 种 拟 周期 波动 。 当 Wi 增 
大 至 3 时 ,波动 幅度 整体 减 小 ,低频 和 次 低频 信号 的 
频率 均 增加 ,因此 波动 幅度 较 大 时 ,频率 较 小 ,波动 
幅度 较 小 时 ,频率 较 大 ; 当 Wi OK 84 以 上 时 ,速度 
波动 消失 ,说 明 流 场 的 不 稳定 只 存在 于 中 等 Wi 数 的 
情况 , 当 弹 性 过 小 或 过 大 时 , 即 完全 黏 性 主导 或 者 完 
全 弹性 主导 时 , 均 不 会 出 现 流动 的 不 稳定 。 

对 胶 束 断裂 指数 X 对 速度 波动 的 影响 进行 了 研 
究 , 如 图 $(c) ~5(d) 所 示 。 发 现 X 减 小 ,速度 波动 
增强 ,波动 信号 的 频率 数 增加 ,并 且 次 级 频率 的 波动 
幅度 增 大 , 即 波动 信号 的 随机 性 增强 ;但 是 当 X 进 一 
步 减 小 时 ,由 流 变 性 质 可 知 ,VCM 模型 逐渐 变 得 与 
UCM 接近 ,难以 获得 收敛 解 , 可 以 推断 速度 波动 随 
着 X 进一步 减 小 呈现 较 强 的 随机 效应 ,这 与 ZHANG 
等 ' 在 实验 中 观察 到 的 颗粒 沉降 速度 变化 特征 定 
性 上 吻合 。 对 于 球 非 稳 态 沉降 问题 ,可 以 推断 是 由 
流体 拉 伸 应 力 对 弹性 变化 响应 的 强 非 线性 效应 
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(b) VCM 流体 圆柱 绕 流速 度 场 随时 间 的 变化 
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(c) VCM 流 体 和 牛顿 流体 绕 流 圆柱 前 后 (7ZOVBD)=0) 的 速度 分 布 


1.0 Xx/D)=1, Y/H)=0 
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(d) 流 场 中 点 (X(x/q)=1 m 的 速度 随时 间 的 变化 
图 4 VOM 流体 圆柱 绕 流 流 场 特征 
(Re «10 ?, BR =0.5 X 20.01,Wi 22) 
Fig.4 The flow field of VCM fluids flowing around a circular 
cylinder (Re ~ 10? , BR 20.5 X 20.01,Wi 22) 
造成 上 述 现 象 的 原因 是 当 Wi de 1. EA PST ds 
黏 性 起 主导 作用 ,其 应 力 分 布 与 在 牛顿 流体 中 分 布 
特征 一 致 (图 6c) ;而 当 Wi » 1 时 ,弹性 作用 增强 ,应 
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力 整体 增 大 , 绕 流 时 圆柱 周围 和 尾 迹 区 剪 切 和 拉 伸 
速率 在 空间 分 布 不 均匀 ,由 Wi = A uy 可 知 绕 流 时 流 
场 不 同位 置 处 的 弹性 大 小 不 同 , 大 应 变 率 区 域 弹 性 
更 强 , 而 由 流 变 曲 线 图 3 可 知 , 拉 伸 应 力 随 着 弹性 大 
小 变化 极 不 均匀 ,导致 应 力 在 空间 上 的 分 布 不 均匀 ， 
随时 间 变 化 趋势 不 一 致 ,产生 流动 不 稳定 。 而 当 Wi 
较 大 时 ,弹性 起 主导 作用 , 避 开 了 应 力 随 弹性 变化 的 
强 非 线 性 区 ,因此 应 力 变 化 比较 平滑 ,不 会 产生 流动 
不 稳定 现象 。 

此 外 研究 了 不 同 的 弹性 和 壁面 约束 下 的 主 应 力 
差 PSD 和 长 链 数量 密度 n, 的 变化 。 主 应 力 差 PSD 
定义 为 Ao = (rs 7 74). +47%, 可 以 描述 胶 束 
在 空间 的 取向 和 伸展 情况 ,与 流动 双 折 射 率 之 差 具 
有 Ac = An/C 的 线性 关系 ,C 为 应 力 光 学 常数 。 图 
6(a) 为 圆柱 绕 流 的 流动 双 折 射 实验 结果 ,图 6 (b) 2g 
Wi 22 时 主 应 力 差 的 模拟 值 ,可 以 看 出 ,模拟 结果 能 
够 较 好 地 反映 出 绕 流 过 程 中 圆柱 周围 胶 束 的 向 列 结 
构 分 布 特性 和 应 力 特 征 ,这 与 SUN 等 2 在 实验 中 观 
测 到 的 圆柱 后 缘 处 的 三 角 区 域 吻合 较 好 。 由 
图 6(f) .图 6(i) 看 出 ,BR 很 小 , 即 壁面 效应 很 弱 时 ， 
在 圆柱 前 缘 处 PSD 最 大 , 尾 迹 区 较 小 ,同样 弹性 (天 
相等 ) 时 ,不 发 生 速度 场 的 不 稳定 ,这 是 因为 尾 迹 区 
应 力 值 大 大 减 小 ,不 会 诱发 流动 的 不 稳定 。 

此 外 通过 对 圆柱 绕 流 后 的 流 场 研究 , 发 现在 雷 
诺 数 Re ~0.01 时 , 绕 流 尾 迹 区 出 现 一 对 上 下 对 称 的 
涡 ( 图 7) ,流体 在 圆柱 表面 发 生 流动 分 离 , 并 且 涡 的 
形态 和 大 小 以 流 场 波动 的 主 频率 作 周 期 性 变化 。 洛 
着 圆柱 后 缘 应 力 值 较 低 的 三 角 区 域内 出 现 一 对 扁平 
的 涡 , 随 着 时 间 增 加 , 涡 沿 着 流动 方向 被 拉 长 ,然后 
又 缩短 恢复 至 初始 的 扁平 状 。 流 体 之 所 以 在 低 雷 诺 
数 绕 流 时 会 从 圆柱 表面 脱落 ,是 因为 在 圆柱 上 下 表 
面 和 尾 迹 区 内 ,流体 剪 切 和 拉 伸 变形 较 大 ,长 胶 束 和 
短 胶 束 均 处 于 伸 长 状态 ,圆柱 表面 薄 层 内 的 向 列 结 
构 的 胶 束 向 下 游 的 过 程 中 为 保证 拉 伸 状态 ,不 会 达 
到 后 驻 点 ,而 是 倾向 于 提前 从 圆柱 表面 脱落 进入 尾 
迹 区 ,进入 尾 迹 区 的 部 分 短 链 由 于 拉 伸 作用 减 小 会 
重新 结合 成 长 链 (图 6g) ,圆柱 后 方 的 三 角 区 域 边 缘 
附近 长 链 数 量 密度 由 明显 增加 , 重 构 后 的 长 链 松 弛 
后 回 缩 ， 重 新 回 到 自由 缠 结 状态 ， 此 区 域内 流体 流 
动 速度 小 ， 部 分 流体 卷 人 到 圆柱 后 方 的 清 止 点 附 
近 ， 形 成 一 对 上 下 对 称 的 涡 ， 因 此 涡 的 形成 是 胶东 
拉 伸 形成 向 列 结 构 及 尾 迹 区 胶 束 的 强 拉 促 共 同 作用 
导致 。 
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5 圆柱 下 游 点 (X=1, 了 =0) 点 速度 变化 的 时 间 谱 和 功率 谱 (Re - 10 ^, BR =0.5) 
Fig.5 The time spectrum and power spectrum of the velocity with the time at point (X 2 1,Y 20) 
behind the cylinder (Re ~ 10? , BR 20.5) 
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0.05 010 0.15 020 023 0.16 0.40 0.60 0.80 1.00 
PSD/Pa 


© 


EL — 84 


(e) 主 应 力 差 PSD 模 拟 值 (Wi=5, BR-0.5) h Jc s ict En 模拟 值 (Wi=5, BR=0.5) 
0 001 002 003 0.044 el 0.94 0.96 0.98 1.00 
PSD/Pa 


O N 


C (M 


(f) 主 应 力 差 PSD 模 拟 值 (Wi=5, BR-0.01) 长 胶 束 数量 密度 n, 模 拟 值 (Wi=5, BR=0.01) 
图 6 绕 流 流体 中 主 应 力 差 (PSD) 和 长 胶 束 数量 密度 分 布 (X=0.01,Re «10 7) 
Fig.6 The principal stress difference and the number density of the long micelle within the fluids (X=0.01,Re~10 s 
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1-10.30 s 
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=. 


E7 绕 流 尾 迹 区 随时 间 变 化 的 对 称 涡 (X=0.005,BR 20.5) 
Fig.7 The symmetrical vortex changing with time in the wake region (X =0.005,BR =0.5) 


牛顿 流体 中 , 当 雷 度数 小 于 0.5 时 , 绕 流 圆柱 
Em 流体 不 会 从 圆柱 壁面 脱落 , 当 雷 诺 数 大 
于 2 时 ,圆柱 后 缘 附近 出 现 一 对 关于 了 = 0 平面 对 称 
的 涡 , 当 雷诺 数 大 于 40 时 , 涡 开 始 变 得 不 稳定 。 而 
在 蠕虫 状 胶 束 流体 中 ,即使 在 很 小 的 雷诺 数 下 ,产生 
了 不 稳定 的 涡 。KHAN 等 呈 也 发 现 了 在 低 雷 诺 数 
(Re ~10“) 和 中 等 维 森 贝 格 数 (天 ~2.5) F, RR 
断裂 因子 上 =0.1 (=0.26) 时 圆柱 后 缘 表面 附近 出 
现 了 一 对 对 称 涡 , 但 是 涡 的 大 小 和 形态 保持 稳定 , 当 
时 ,并 未 出 现 流 场 的 非 稳 态 波动 ,因此 涡 的 

态 和 大 小 也 不 变 , 当 & 减 小 至 0.01 时 ， 出 现 流 声 
5. 涡 的 形态 和 大 小 也 随 之 发 生变 化 ,这 
与 本 研究 模拟 结果 基本 一 致 。 上 表征 长 胶 束 链 断 裂 
的 非 线 性 效应 大 小 ,é 越 小 , 胶 束 拉 伸 应 力 对 弹性 响 
应 的 非 线 性 效应 越 显著 , 越 容易 出 现 流 动 不 稳 定 , 流 
场 的 不 稳定 导致 涡 的 形态 和 大 小 随时 间 变 化 ,é 越 
大 , 胶 束 拉 伸 应 力 对 弹性 响应 的 非 线 性 效应 越 弱 , 越 
不 容易 发 生 流动 不 稳定 , 当 & 较 大 时 ,流动 呈 稳 态 ， 
涡 的 大 小 和 形态 也 不 随时 间 变 化 。 

此 外 文献 中 采用 Oldroyd-B 和 Giesekus 等 本 构 
模型 时 ,并 未 发 现在 低 雷 诺 数 时 的 流 场 不 稳定 现 
45/755] mi OLIVEIRA 4&^* 3x JH FENE-CR 本 构 模 
型 研究 圆柱 绕 流 时 发 现 ,完全 忽略 惯性 效应 ， 
大 于 1.3 时 ,圆柱 后 缘 发 现 了 形态 和 大 小 周期 性 
化 的 涡 , 而 本 研究 得 到 的 流 场 波动 具有 拟 周期 性 3 
且 当 拉 伸 应 力 随 弹性 变化 的 非 线 性 较 强 时 ,波动 非 
周期 性 增强 ,采用 与 OLIVEIRA 等 ”” 同样 的 流 变 参 
数 求解 其 流 变 曲线 ,发 现 Wi 从 0.1s EI ls, 
拉 伸 黏度 增 大 了 50 倍 以 上 , 拉 伸 应 力 随 弹性 出 现 较 
强 的 非 线性 变化 趋势 ,而 剪 切 应 力 随 弹性 呈 线 性 变 
化 ,因此 进一步 说 明 拉 伸 应 力 随 弹性 变化 的 强 非 线 
性 效应 导致 流动 不 稳定 现象 的 发 生 ,并 且 由 于 足 


够 小 时 VCM 模型 拉 伸 应 力 随 弹性 变化 的 非 线 性 比 
FENE-CR 模型 更 强 , 流 场 波 动 的 幅度 较 大 ,并 且 波 
动 非 周期 性 更 强 。 


3 结 论 


采用 双 组 分 VOM 黏 弹性 流体 本 构 方 程 描述 蠕 
虫 状 胶东 溶液 的 黏 弹性 和 胶 束 链 断 裂 ,研究 了 其 绕 
圆柱 流动 时 流动 状态 以 及 流 变 参数 对 流动 不 稳定 性 
的 影响 。 结 论 如 下 。 

1) 低 雷诺 数 蠕 虫 状 胶 柬 涂 液 圆柱 绕 流 中 ,在 中 
等 维 森 贝 格 数 下 可 能 发 生 流 动 不 稳定 ,流动 完全 由 
弹性 主导 或 者 黏 性 主导 时 均 不 会 出 现 流动 不 稳定 
现象 。 

2) 在 低 雷 诺 数 下 , 圆柱 后 缘 发 生 流动 分 离 , 出 
现 一 对 周期 性 变化 的 对 称 涡 ,说 明 流体 弹性 会 促进 
流动 分 离 和 流动 不 稳定 的 发 生 ,流动 分 离 是 由 黏 弹 
性 胶 束 链 的 向 列 分 布 结构 和 尾 迹 区 的 强 拉 伸 共 同 作 
用 引起 。 

3) 流 体 拉 伸 应 力 对 弹性 的 强 非 线 性 响应 是 导致 
流 场 不 稳定 的 根本 原因 , 拉 伸 应 力 随 弹性 变化 的 非 
线性 效应 较 弱 时 ,速度 场 呈 稳定 状态 , 拉 伸 应 力 随 弹 
性 变化 的 非 线 性 较 强 时 ,出现 流 场 的 不 稳定 , 拉 伸 应 
力 随 弹 性 变化 的 非 线性 越 强 , 速度 场 波 动 不 规则 性 
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